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The  turbulent  f lowf  ield  in  a backmixed  combustor  is  modeled  analytically  by  simultaneously 
solving  the  governing  partial  differential  conservation  equations.  An  algebraic  and  a two-equa- 
tion eddy  viscosity  model  are  employed  in  the  numerical  simulation  to  account  for  the 
turbulence  transport  processes.  Predicted  distributions  of  isothermal  flow  properties  are 
systematically  compared  to  experimentally  obtained  data  to  appraise  the  eddy  viscosity  models. 
In  addition,  predicted  distributions  of  reacting  flow  properties  are  presented  to  illustrate 
the  applicability  of  current  numerical  methods  to  the  prediction  of  continuous  combustion 
flows. 

The  turbulent  momentum  and  mass  transport  properties  are  evaluated  in  isothermal  flow 
for  a range  of  mixing  conditions.  Both  eddy  viscosity  models  qualitatively  describe  the  system 
hydrodynamics,  but  the  detailed  flow  structure  is  inadequately  represented.  The  mass  transport 
predictions  from  the  algebraic  viscosity  model  agree  favorably  with  experiment.  The  inferior 
performance  of  the  two-equation  viscosity  model  is  remedied  by  refining  the  boundary 
condition  specification  for  the  turbulence  energy  dissipation  rate.  It  is  shown  that  the  turbulence 
energy  dissipation  rate  adjacent  to  critical  solid  walls  strongly  influences  the  overall  mixing 
characteristics  of  the  two-equation  model. 

The  isothermal  flow  results  indicate  that  the  algebraic  eddy  viscosity  model  provides 
cost-effective  predictions  of  the  general  fluid  flow  patterns  and  mass  transport  trends  in 
confined  flows  exhibiting  strong  recirculation.  The  two-equation  eddy  viscosity  model  provides 
better  resolution  of  the  small-scale  turbulence  processes  but  requires  careful  testing  to  ensure 
realistic  predictions. 

The  hot  flow  calculations  accentuate  the  inadequate  transport  characteristics  identified 
in  the  isothermal  flow  analysis  and  verify  that  considerable  testing  of  the  numerical  model 
is  required  before  proceeding  to  the  complicating  conditions  of  combustion. 


Introduction 


The  design  and  operation  of  continuous 
combustion  devices  such  as  gas  turbines, 
boilers,  and  furnaces  may  be  assisted  by  the 
development  and  application  of  suitable  pre- 
dictive models  that  account  for  pollutant  pro- 
duction, combustion  efficiency,  and  heat  re- 
lease behavior.  The  emergence  of  continuum 
flow  models  formulated  from  governing  con- 
servation equations  represents  a significant 
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advance  toward  the  numerical  prediction  of 
combustipn-chamber  pierformance.  The  ap- 
plication of  continuum  models  to  the  design 
of  combustion  hardware  is  attractive  because 
the  point-by-point  prediction  of  velocity,  tem- 
perature, and  composition  is  an  integral  part 
of  the  computational  procedure.  The  quantita- 
tive prediction  of  turbulent,  backmixed  c-om- 
bustion  has  not  yet  been  achieved  because  of 
uncertainties  linked  to  the  submodels  of 
fundamental  combustion  processes.  The  pres- 
ent study  addresses  an  essential  element  of 
continuum  flow  modeling — the  use  of  eddy 
viscosity  submodels  to  account  for  turbulence 
transport  processes.  The  objective  is  to  deter- 
mine the  utility  and  applicability  of  eddy 
viscosity  submodels  in  the  numerical  simula- 


1675 


i . 

. t 

1676  TURBULCNT  flames  AND  COMBUSTION 


tion  of  flows  exhibiting  strong  recirculation 
or  backmixing. 


Approach 

Selected  eddy  viscosity  submodels  are 
systematically  evaluated  in  a combined  ana- 
lytical and  experimental  investigation.  A back- 
mixed  laboratory  combustor  is  employed  to 
conduct  tests  for  both  isothermal  (nonreacting) 
and  hot  (reacting)  flow  conditions.  For  the 
case  of  isothermal  flow,  the  eddy  viscosity 
submodels  are  first  evaluated  for  the  prediction 
of  system  hydrodynamics.  A tracer  gas  is 
subsequently  introduced  to  evaluate  the  eddy 
viscosity  submodels  for  nonuniform,  multi- 
component  flow  where  mass  as  well  as  mo- 
mentum exchange  is  important.  Experimental 
measurements  and  numerical  predictions  of 
velocity  and  tracer  isopleths  for  a range  of 
approach  velocities  and  backmix  charac- 
teristics provide  the  critical  tests. 

The  results  reported  here  emphasize  the  case 
of  isothermal  flow.  A satisfactory  assessment 
and  understanding  of  the  isothermal  case  is 
a prerequisite  to  engaging  the  hot  flow  case. 
Combustion  imposes  the  additional  require- 
ment for  submodels  of  chemistry  and  turbu- 
lence/chemistry interaction  that  may  mask 


fundamental  deficiencies  in  the  eddy  viscosity 
submodel.  However,  hot  flow  results  are  pre- 
sented as  an  indication  of  the  suitability  of 
the  two  eddy  viscosity  submodels  in  the  pre- 
diction of  reacting  flows. 

The  results  provide  useful  information  for 
the  development  of  applied  computational 
tools  for  analyzing  the  performance  of  contin- 
uous combustion  systems.  In  particular,  the 
results  identify  strengths  and  weaknesses  in 
the  use  of  two  eddy  viscosity  submodels,  and 
demonstrate  the  need  to  thoroughly  test  pre- 
dictive methods  for  the  case  of  isothermal  flow 
before  proceeding  to  the  case  of  hot  flow. 


Combustion  System 

The  experimental  configuration  is  a 51  mm 
I.D.  X 457  mm  cylindrical  Vycor  combustion 
chamber  containing  an  aerodynamic  (oppos- 
ed-jet) flameholder  as  shown  in  Fig.  1.  The 
incoming  flow  of  premixed  methane  and  air 
is  opposed  by  a high  velocity  jet  (.rh./rh^  <<  1) 
issuing  from  a 1.3  mm  I.D.  (^4  mm  O.D.) 
water-cooled  tube.  The  jet  is  coincident  with 
the  combustor  axis  and  located  80  mm  up- 
stream from  the  combustor  exit.  The  back- 
mixed  zone  required  for  flame  stabilization 
is  generated  along  the  jet  boundary. 


-FLAME  FRONT  pRECIRCULATING,  PARTIALLY 
OR  COMPLETELY  BURNT  GAS  , 


STAGNATION 

POINT 


Fic.  1.  Opposed-jet  combustor  flowfield  schematic 
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The  opposed-jet  combustor  (OJC)  exhibits 
essential  features  found  in  practical  continu- 
ous combustion  systems  (backmixing,  high 
intensity  combustion)  but  avoids  the  attendant 
complexities  (nonpremixed,  two  phase  com- 
bustion) appropriate  for  subsequent  studies. 
The  separation  of  the  critical  backmixing  re- 
gime from  influential  boundary  conditions  is 
important  to  testing  the  eddy  viscosity  sub- 
models. In  addition,  the  OJC  provides  a ver- 
satile flowfield  to  challenge  the  predictive 
method.  For  example,  the  rate  and  intensity 
of  backmixing  in  the  recirculation  zone  may 
be  altered  experimentally  and  numerically  over 
a broad  range  by  varying  the  ratio  of  approach 
velocity  to  jet  velocity. 


This  approach  is  standard  practice  in  engi- 
neering calculations  of  turbulent  flows  and 
encompasses  the  introduction  of  turbulent  ex- 
change coefficients  that  imply  Newtonian 
stress-strain  relationships,  Fourier  heat  con- 
duction, and  Fickian  diffusion.  A turbulent 
“eddy”  viscosity,  p.,,  is  used  to  prescribe  the 
turbulent  momentum  flux,  i.e. 


The  turbulent  transport  of  scalar  quantities 
is  related  to  the  eddy  viscosity  through  the 
appropriate  Prandtl  or  Schmidt  numbers  as 
follows: 


Numerical  Method 

The  theoretical  analysis  of  the  opposed-jet 
combustor  is  accomplished  by  solving  nu- 
merically the  time-averaged  conservation 
equations  for  turbulent,  backmixed  flow  with 
provision  for  chemical  reaction.  The  backmix- 
ing requires  that  the  governing  mass,  momen- 
tum, energy,  and  species  conservation  equa- 
tions be  elliptic.  The  computational  procedure 
used  to  solve  simultaneously  the  elliptic  partial 
differential  equations  for  the  backmixed  flow- 
field  is  based  upon  the  TEACH  method  de- 
scribed by  Gosman  et  al.‘  Equations  contain- 
ing the  primitive  (p  - V)  variables  are  solved 
simultaneously  to  generate  an  isothermal  flow 
solution.  Equations  for  stagnation  enthalpy  (ft) 
and  species  mass  fraction  (Y,)  are  then  added 
to  obtain  a solution  for  hot  flow  conditions. 
The  dependent  variables-  (<|>)  are  cast  into  a 
generalized  equation  of  the  form: 


The  present  study  explores  specification  of 
the  eddy  viscosity,  p,,  and  turbulent  mass 
transport,  in  confined  flows  exhibiting 
strong  backmixing  or  recirculation.  Compara- 
tive calculations  are  performed  using  two  dis- 
tirtct  turbulence  submodels  to  define  the  eddy 
viscosity — an  algebraic  model  and  a more 
sophisticated  two-equation  submodel: 

Algebraic* 

p,  = ED*^*L-‘'3p*/3(mi;»)>/»  (4) 

K = empirical  coefficient 
O = diameter  of  combustor 
L = length  of  combustor 
p = local  density 
mU*  ^ inlet  kinetic  energy 


d Id  d / d4t\ 

— (pUb)  + r(pV<j))  = — r,.*— ) 

dx  r dr  dx  \ * dx  / 

1 a / a<|>\ 

r dr  \ dr  / * 

where  r,_^  and  represent  the  turbulent 
exchange  coefficient  and  source  terms  respec- 
tively. Boundary  conditions  are  prescribed  for 
all  dependent  vari^jbles  following  the  recom- 
mendations of  Ref.  (1).  After  specification  of 
inlet  and  solid  boundary  conditions,  the  finite 
difference  calculations  are  carried  out  line-by- 
line over  an  18  X 31  nonuniform  grid  covering 
the  half-plane  of  the  axisymmetric  flowfield. 

Turbulence  Submodels 

Turbulent  transport  properties  are  evaluated 
using  time-mean  or  "effective”  quantities. 


Ttoo-Equation 

P,  *c,,pkV«  (5) 

c^  = empirical  coefficient 
p = local  density 
k = local  turbulent  kinetic  energy 
c = local  turbulent  energy  dissipation 
rate 


The  algebraic  formulation  is  commonly  used 
in  predictive  modeling.  The  variation  in  eddy 
viscosity  depends  only  upon  local  density 
(spatially  invariant  for  isothermal  conditions) 
for  a given  geometry  and  inlet  flow  rate.  The 
two-equation  subm^el  introduces  additional 
dependent  variables  that  complicate  the  nu- 
merics and  add  computational  cost.  However, 
the  two-equation  submodel  accounts  for  spa- 
tial variations  in  turbulent  kinetic  energy  and 
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length  scale,  and  has  been  tested  satisfactorily 
in  a number  of  fluid  dynamic  encounters.^ 


Results  (Isothermal  Flow) 


I 

\ 


To  effectively  test  the  submodels  of  eddy 
viscosity  over  a range  of  turbulent  conditions, 
experimental  data  were  obtained  for  approach 
velocities  of  15.24  and  7.62  m/sec  and  jet 
velocities  of  30.5,  61  and  130  m/sec.  Velocity 
profiles  were  mapped  throughout  the  f lowf  ield 
using  a 1 mm  O.D.  pitot  tube.  Tracer  profiles 
were  generated  by  introducing  alternatively 
pure  carbon  monoxide  (CO)  in  the  jet  and  pure 
carbon  dioxide  (COj)  in  the  mainstream.  The 
mixing  between  the  jet  and  mainstream  was 
assessed  by  mapping  CO  and  COj  profiles 
throughout  the  flowfield  using  a 3 mm  O.D. 
stainless  steel  tube  and  Beckman  315B  NDIR 
instruments  for  the  gas  sampling  and  analysis. 
The  probe  dimensions  are  small  (probe 
O.D./OJC  l.D.  = .06)  to  minimize  hydrody- 
namic interferences.  The  Reynolds  number  of 
the  approach  stream  ranged  from  2.5  x lO'* 
(at  7.62  m/sec)  to  5.0  x 10'*  (at  15.24  m/sec). 


Momentum  Transport 

The  eddy  viscosity  submodels  were  first 
tested  for  a homogeneous,  isothermal  flowfield 
to  explore  the  effectiveness  of  the  submodels 
in  predicting  momentum  transport.  Experi- 
mental and  predicted  profiles  of  the  axial 
component  of  mean  velocity  are  presented  in 
Fig.  2 for  an  average  inlet  velocity  of  15.24 
m/sec.  The  selected  profiles  emphasize  the 
flow  regime  dominated  by  recirculation.  The 
velocity  data  are  subject  to  experimental  errors 
on  the  order  of  five  percent  in  the  bulk  flow 
region,  but  in  the  recirculation  zone  the  results 
are  compromised  by  pitot  tube  errors  of 
twenty-five  percent  or  more.  The  calculated 
velocity  profiles  for  both  turbulence  submod- 
els assume  an  experimentally  determined  ve- 
locity profile  as  an  inlet  condition. 

The  trend  of  the  velocity  profiles  agree 
qualitatively  with  the  experimental  results. 
The  two-equation  submodel  retains  the  inlet 
turbulent  velocity  profile  in  the  approach  sec- 
tion while  the  velocity  profiles  predicted  by 
the  algebraic  submodel  flatten  out  due  to  the 
spatially  uniform  viscosity.  The  presence  of 
the  opposing-jet  is  communicated  further  up- 
stream for  the  two-equation  submodel  than 
for  the  algebraic  eddy  viscosity,  and  the  two- 
equation  submodel  demonstrates  better  cor- 
relation with  the  experimentally  observed 
stagnation  point. 


The  mean  velocity  profiles  in  the  recircula- 
tion zone  are  also  better  represented  by  the 
two-equation  turbulence  submodel,  although 
both  submodels  over-predict  the  jet  expansion. 
In  addition,  the  two-equation  velocity  profiles 
decrease  too  rapidly  near  the  chamber  wall. 
Downstream  from  the  backmixed  zone  both 
submodels  fail  to  adjust  for  the  redistribution 
of  the  bulk  flow  in  the  turbulent  wake  of  the 
recirculation  zone  and  under-predict  the  ve- 
locity adjacent  to  the  jet  wall.  Both  submodels 
display  good  agreement  in  the  bulk  flow  re- 
gime in  the  absence  of  recirculation  or  wall 
effects. 

The  recirculation  zone  strength  was  exam- 
ined by  reducing  the  approach  velocity  by  a 
factor  of  two  while  retaining  the  same  jet 
velocity.  The  experimental  and  predicted  ve- 
locity profiles  are  presented  in  Fig.  3.  The 
calculated  results  demonstrate  a significant 
departure  from  exjjeriment  in  the  region  of 
strong  backmixing.  The  reverse  flow  region 
again  penetrates  further  upstream  in  the  two- 
equation  submodel  prediction. 

The  acceleration  of  the  bulk  flow  around 
the  recirculation  zone  is  under-predicted  in 
contrast  to  the  over-predicted  jet  expansion  for 
both  submodels.  This  result  is  especially  evi- 
dent near  the  plane  of  the  jet  exit  where  the 
redistribution  of  the  bulk  flow  with  the  jet 
discharge  is  restricted.  This  effect  is  conveyed 
downstream  from  the  recirculation  zone  where 
the  two-equation  velocity  profiles  correlate 
well  with  experiment  near  the  chamber  wall 
but  deviate  near  the  jet  wall. 

Mass  Transport 

The  two  eddy  viscosity  submodels  were 
secondly  tested  in  conjunction  with  the  turbu- 
lent mass  exchange  submodel  (Eq.  3)  for  the 
conditions  of  an  isothermal,  nonhomogeneous 
flowfield.  Spatial  distributions  of  the  predict- 
ed tracer  concentrations  are  compared  to  ex- 
perimental measurements  in  Figs.  4 and  5 for 
the  baseline  conditions  outlined  previously. 
The  tracer  gas  concentration  data  were  ob- 
tained at  radial  increments  of  1.3  mm  and  axial 
stations  every  13  mm.  The  concentration  iso- 
pleths  (experimental)  were  constructed  from 
values  measured  at  the  discrete  data  points 
in  the  experimental  flowfield  and  compared 
to  concentration  isopleths  (predicted)  con- 
structed from  values  predicted  at  the  discrete 
grid  points  in  the  computational  field. 

Figure  4 presents  the  algebraic  and  two- 
equation  submodel  predictions,  and  experi- 
mental tracer  conc«ntration  distributions  for 
a 15.2  m/sec  average  inlet  velocity.  (The  jet 
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Fic.  4.  Isothermal  flow  tracer  concentration 
profiles.  0„  = 15.24  m/sec;  0^  = 130m/sec;  <I>„  = 
0.;  100%  CO  Jet. 


stream  is  100  percent  carbon  monoxide.)  Im- 
mediately evident  is  the  contrasting  axial  and 
radial  spread  of  the  predicted  tracer  concentra- 
tion profiles.  Differences  in  axial  transport  of 
CO  are  consistent  with  the  observations  for- 
warded in  the  previous  section  on  momentum 
transport.  More  notable  are  the  diverse  radial 
transport  predictions.  The  spreading  of  the  jet 
discharge  into  the  mainstream  for  the  algebraic 
eddy  viscosity  conforms  to  the  experimental 
results.  In  contrast,  the  two-equation  submodel 
limits  the  radial  mixing  of  the  jet  stream  (tracer) 
with  the  bulk  flow.  The  segregation  of  the 
jet  and  the  mainstream  is  verified  by  the  high 
CO  concentrations  along  the  jet  wall,  in  con- 
trast to  the  absence  of  CO  near  the  chamber 
wall.  This  result  was  further  substantiated  by 
introducing  a carbon  dioxide  (CO^)  tracer  in 
the  approach  stream.  The  numeric^  analysis 
of  CO.  concentrations  for  the  two-equation 
submodel  predict  negligible  mass  transport  of 
CO,  from  the  mainstream  into  the  recirculation 
zone.  Downstream,  the  predicted  mixing  of 
CO,  from  the  bulk  stream  into  the  flow  asso- 
ciated with  the  jet  stream  is  minor.  Experi- 
mentally CO,  was  observed  to  be  well  mixed 
into  the  jet. 


a)  M'OlQCbraic 


biM'two -equation 


Fic.  5.  Isothermal  flow  tracer  concentration 
profiles.  = 7.62  m/sec;  Oj  = 130  m/sec;  <I>„ 
= 0.;  100%  CO  Jet. 


The  results  for  the  strong  recirculation  zone 
(Fig.  5)  demonstrate  similar  trends.  Although 
the  velocity  data  presented  earlier  in  Fig.  3 
gave  poorer  agreement  for  conditions  of  strong 
recirculation,  the  mixing  is  improved  because 
of  the  enhanced  backmixing. 


Results  (Hot  Flow) 

A two-step  global  reaction  mechanism  for 
methane  oxidation  was  adopted  for  the  present 
study  to  illustrate  the  applicability  of  current 
numerical  methods  to  the  prediction  of  contin- 
uous combustion  flows,  and  to  briefly  explore 
the  suitability  of  the  eddy  viscosity  subm^els 
for  the  case  of  hot  flow: 

CH^ -I- 0,-»CO -H  H,0  (6) 

CO  -I-  O,  -►  CO,  (7) 

A detailed  description  of  the  numerical  for- 
mulation of  the  reaction  rate  expressions  and 
boundary  conditions  may  be  found  elsewhere.* 
The  dependent  variables  for  the  hydrocarbon 
system  include  the  mass  fractions  and 

^co*-  Distributions  of  other  major  species — 
water  (H,0),  oxygen  (O,)  carbon  monoxide 
(CO) — are  related  to  ^co, 

mental  mass  conservation.  An  initial  solution 


1682 


n KBULENT  Fl^MES  A.\0  COMBUSTION 


for  the  reacting  fiowfielcl  was  oirtained  using 
published  reaction  rates'  which  were  later 
refined  to  more  closely  resemble  the  experi- 
menta  data.  This  approach  was  justified  on 
the  basis  that  the  reaction  rate  data  were 
obtained  from  a well-stirred  reactor  study  and 
may  not  apply  to  the  cm  rent  study,  which 
deviates  substantially  from  well-stirred  condi- 
tions. Global  reaction  rates  were  judged  ac- 
ceptable for  the  present  study  Ixjcause  of  the 
uncertainties  encompassing  the  specification 
of  the  system  aerodynamics. 

Since  the  aim  of  the  study  is  to  predict  the 
performance  of  backmixed  combustion  sys- 
tems, the  hot  flow  calculations  were  extended 
to  include  the  prediction  of  the  pollutant 
species,  nitric  oxide  (NO).  Nitric  oxide  kinetics 
were  based  on  the  familiar  Zeldovich'*  mecha- 
nism: 

O-1-N2— NO-l-N  (8) 

N -h  Oj  ->  iNO  -H  O (9) 

Noting  that  reaction  (8)  is  the  rate  limiting 
step  and  adopting  the  simplifying  assumption 
of  O/Oj  etjuilibrium  results  in  the  reaction 
rate  expression: 


t/[NO] 

^ = 2fc„K„[NJ[Oj‘'"  (10) 

at 

The  temperature  and  Oj  distributions  obtained 
from  the  solution  of  the  hydrocarbon  systems 
were  used  as  a basis  for  the  NO  kinetic  cal- 
culations. 

The  experimental  tests  conducted  to  com- 
plement the  numerical  predictions  of  hot  flow 
properties  utilized  stoichiometric  proportions 
of  premixed  methane  and  air  in  Ixrth  the  main 
and  jet  streams.  The  reactants  were  initiall>' 
at  ambient  temperature  and  the  combustion 
was  completed  at  atmospheric  pressure.  The 
test  matrix  included  the  same  approach  veloci- 
ties selected  for  the  isothermal  cases — 7.62  and 
15.24  m/sec — and  a jet  velocity  of  130  m/sec. 

Selected  results  of  hot  flow  properties  are 
presented  in  Figs.  6,  7,  and  8 for  the  conditions 
indicated.  Velocity,  temperature,  and  NO  data 
are  presentetl  as  indicators  of  useful  design 
information  that  may  be  derived  from  the 
numerical  simulation  of  backmixed  combus- 
tion processes.  A phrrtograph  of  the  opposed-jet 
combustor  is  shown  in  Fig.  6 with  predicted 
velocity  vectors  superimposed.  The  detailed 
maps  of  ‘^lowfield  properties  pre.sented  in  Figs. 


Fit;.  6.  Photograph  of  opposed-jet  combustor  burning  premixed  methane  and  air  with  predictetl  velocity 
vectors  siiperimposetl.  0„  = 7.62  m/sec;  0^  = 1.30  ni/sec;  <l>„  = 4>j  = 1. 
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Fic.  7.  Hot  flow  temperature  distribution.  (7„  = 
15.24  m/sec;  Oj  = 130  m/sec;  4>„  = = 1.0. 


b)  M'  two-equQlion 


Fig.  8.  Hot  flow  nitric  oxide  concentration  pro- 
files. « 15.24  m/sec;  Of  “ 130  m/sec; 
«<!>(«  1.0. 


7 and  8 were  obtained  by  conventional  probing 
techniques.  Temperature  measurements  were 
made  using  an  uncoated  Pt/Pt-13%  Rh  fine- 
wire  thermocouple.  The  data  are  not  corrected 
for  radiation  losses,  but  provide  an  adequate 
description  of  the  OJC  heat  release  distribu- 
tion. A Scott  Model  125  chemiluminescent 
analyzer  was  used  for  oxides  of  nitrogen  mea- 
surements. Gas  samples  were  extracted  via  a 
moveable  3.2  mm  O.D.  watercooled,  316  stain- 
less steel  probe  having  a tubular  inlet,  and 
conveyed  through  a heated  teflon  sample  line 
to  a packaged  exhaust  gas  analysis  system.  Hie 
measurement  of  species  susceptible  to  sample 
transformations  (e.g.  NO,  NOj)  must  be  re- 
garded as  qualitative  estimates  of  the  actual 
local  concentrations. 


Discussion 

Isothermal  Flow 

The  results  of  the  momentum  and  mass 
transport  studies  identify  several  important 
characteristics  of  the  turbulence  submodels. 
The  turbulent  viscosity  from  the  two-equation 
submodel  varies  through  the  flowfield  as  is 
physically  expected.  The  two-equation  sub- 
model also  demonstrates  the  generation  and 
dissipation  of  turbulence  kinetic  energy  near 
the  jet  exit,  the  stagnation  point  and  near  the 
boundaries  consistent  with  other  numerical 
investigations.^  In  contrast,  the  turbulent  vis- 
cosity from  the  algebraic  submodel  is  spatially 
uniform  throughout  the  flowfield  for  isother- 
mal flow,  and  varies  only  with  a change  in 
the  kinetic  energy  of  the  inlet  streams.  The 
impact  of  these  differences  between  the  two 
eddy  viscosity  submodels  is  evidenced  by  the 
velocity  profiles,  the  predicted  location  of  the 
stagnation  point,  the  radial  mass  transport,  and 
the  dissipation  of  the  kinetic  energy  of  turbu- 
lence in  the  wake  region. 

Velocity  Profiles.  The  predicted  velocity 
profiles  agree  qualitatively  with  experimental 
results.  The  two-equation  submodel  more 
closely  approximates  the  experimental  trends 
although  the  predicted  velocity  contours  for 
conditions  of  strong  recirculation  are  quantita- 
tively incorrect  throughout  the  flowfield  for 
both  turbulence  submodels.  The  departure 
from  experiment  is  especially  evident  near 
walls,  and  suggests  the  need  to  refine  the  wall 
functions  used  to  specify  boundary  conditions 
for  turbi  (“nt  kinetic  energy  or  dissipation  rate, 
and/or  to  decrease  the  grid  spacing  near  the 
wall. 

Stogruition  Point.  The  location  of  the 
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stagnation  point  as  indicated  by  the  velocity 
profiles  and  upstream  extent  of  tracer  mea- 
surements is  better  described  by  the  two-equa- 
tion submodel.  The  effective  viscosity  predict- 
ed by  the  algebraic  submodel  is  spatially 
uniform  through  the  flowfield  for  given  inlet 
conditions.  No  account  is  made  for  changes 
in  the  kinetic  energy  of  turbulence  or  energy 
dissipation  rate  in  the  region  of  encounter 
between  the  main  and  jet  flows.  The  lower 
viscosity  predicted  in  the  recirculation  zone 
by  the  two-equation  submodel  is  more  effec- 
tive in  describing  the  actual  axial  momentum 
and  mass  exchange. 

Radial  Mass  Transport.  Although  the  ex- 
perimentally determined  tracer  isopleths  are 
in  closest  agreement  in  the  recirculation  zone 
for  the  two-equation  submodel,  the  down- 
stream radial  spread  is  more  effectively  de- 
scribed by  the  algebraic  submodel.  The  uni- 
formly high  viscosity  predicted  by  the  alge- 
braic turbulence  submodel  aids  the  overall 
mass  transport  throughout  the  flowfield,  i.e. 
axial  and  radial  diffusion  from  the  jet  is  greater 
for  the  algebraic  rather  than  for  the  two-equa- 
tion submodel. 

The  radial  mass  transfer  predicted  by  the 
two-equation  turbulence  submodel  suggests 
that  the  principal  mechanism  of  mass  transport 
is  by  large-scale  convection  rather  than  by 
small  scale  gradient  diffusion  processes.*”  The 
tracer  concentration  varies  gradually  through 
the  recirculation  zone.  The  downstream  mix- 
ing layer  is  bounded  by  steep  concentration 
gradients  that  enclose  the  jet  mixing  region. 
The  numerical  predictions  indicate  that  at  high 
approach  velocities  the  bulk  flow  detours 
around  the  recirculation  zone,  retaining  up- 
stream flow  properties  (low  viscosity  and 
small-scale  mixing)  and  effectively  precludes 
radially-directed  mass  transport.  Conversely, 
the  experimental  trend  reinforces  the  conjec- 
ture that  turbulent  exchange  processes  should 
be  active  throughout  the  annulus  downstream 
from  the  recirculation  zone. 

Wake  Region.  The  kinetic  energy  of  turbu- 
lence predicted  by  the  two-equation  submodel 
is  quickly  dissipated  in  the  wake  region  down- 
stream of  the  recirculation  zone.  The  resulting 
low  viscosity  impedes  turbulent  mixing 
throughout  the  wake  region.  In  contrast,  the 
uniformly  high  viscosity  of  the  algebraic 
submodel  is  effective  in  predicting  the  radial 
spread  of  tracer  observed  experimentally. 

Dioffiosis 

The  results  and  subsequent  discussion  have 
exposed  deficiencies  in  both  of  the  eddy  vis- 


L 


cosity  submodels  that  were  tested  and  have 
identified  directions  for  possible  refinements. 
A supplementary  investigation  was  initiated 
to  determine  the  origin  of  the  shortcomings 
of  the  two-equation  eddy  viscosity  submodel 
in  the  present  application  and  to  improve  the 
IJerformance  of  the  submodel  in  strongly  re- 
circulating flows. 

A series  of  numerical  experiments  were  con- 
ducted to  examine  the  effect  of  boundary 
condition  specifications  on  the  turbulence 
transport  characteristics  of  the  two-equation 
eddy  viscosity  submodel  for  isothermal  flow 
conditions.  Turbulence  properties  at  the  jet 
and  main  stream  inlets  and  adjacent  to  solid 
walls  were  systematically  varied  to  improve 
numerical  /experimental  correlation  and  to  de- 
termine the  sensitivity  of  the  solution  to  inlet 
and  solid  boundary  conditions. 

Modifications  to  the  main  inlet  conditions, 
i.e.  increased  turbulence  kinetic  energy  ( k ) and 
decreased  turbulence  energy  dissipation  (e), 
improved  the  radial  spread  of  tracer  con- 
centration but  did  not  overcome  the  experi- 
mental discrepancies.  The  increased  viscosity 
in  the  main  stream  also  tended  to  move  the 
predicted  stagnation  point  downstream.  Simi- 
lar variations  at  the  jet  inlet  had  little  impact 
on  the  radial  mass  transport  but  had  a signifi- 
cant influence  on  the  predicted  stagnation 
point.  Parametric  changes  reducing  the  viscos- 
ity in  the  vicinity  of  the  jet  exit  tended  to 
extend  the  stagnation  point  upstream. 

The  principal  outcome  of  the  parametric 
studies  was  the  observation  that  in  confined 
flows,  where  the  recirculation  zone  may  c-on- 
tact  solid  boundaries,  the  specification  of 
near-wall  turbulence  energy  dissipation  rate 
may  significantly  influence  the  mixing  charac- 
teristics of  the  two-equation  turbulence 
submodel.  These  results  are  consistent  with 
the  findings  of  earlier  studies  where  a reduc- 
tion of  the  energy  dissipation  rate  adjacent 
to  a downstream  facing  wall  benefited  nu- 
merical/experimental correlation.  Similarly, 
the  reduction  of  the  energy  dissipation  rate 
on  the  upstream  facing  step  of  the  jet  tube 
was  found  to  radically  alter  the  predicted  mass 
transport  behavior  of  the  two-equation  sub- 
model. The  improved  radial  mass  transport 
is  evidenced  by  the  tracer  concentration  profi- 
les in  Fig.  9.  The  stagnation  point  correlation 
may  be  restored  by  modifying  the  turbulence 
kinetic  energy  or  tbe  rate  of  turbulence  energy 
dissipation  at  the  jet  inlet. 

The  analysis  has  shown  that  the  turbulence 
energy  dissipation  rate  in  the  flowfield  is 
highly  dependent  on  the  shear  stress  along 
critical  solid  walls.  The  initially  p>oor  perform- 
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Fic.  9.  Isothennal  flow  tracer  concentration 
profiles  for  refined  jet  boundary  conditions  (p.: 
two-equation).  0„  = 15.24  m/sec;  t)j=  130  m /sec; 
<b„  = 0.;  100%  CO  Jet. 

ance  of  the  two-equation  eddy  viscosity 
submodel  was  remedied  by  refining  the  wall 
boundary  conditions  for  the  turbulence  energy 
dissipation  rate.  The  supplementary  investiga- 
tion has  also  demonstrated  the  ne 'd  to  test 
the  two-equation  turbulence  submodel  for 
isothermal  flow  conditions  prior  to  simulating 
hot  flows. 

Hot  Flow 

The  predicted  hot  flow  results  for  the  base- 
line eddy  viscosity  submodels  qualitatively 
describe  the  bulk  convective  flowfield  and  the 
flameholding  characteristics  of  the  OJC,  but 
the  extent  of  chemical  reaction  is  generally 
overestimated.  The  accelerated  consumption 
of  fuel  and  oxygen  and  the  assumed  adiabatic 
boundary  conditions  precipitate  elevated 
flame  temperatures  and  excess  nitric  oxide 
production.  The  spatial  distributions  of  hot 
flow  properties  substantiate  the  inherent  limi- 
tations of  the  eddy  viscosity  submodels  that 
emerged  in  the  isothermal  flow  analysis.  The 
poor  numerical  simulation  of  the  detailed 
flame  structure  is  attributed  to  deficiencies  in 
the  coupled  turbulence/chemistry  submodels 
and  the  boundary  condition  specifications. 


Conclusions 

The  relative  merits  of  an  algebraic  and  a 
two-equation  eddy  viscosity  submodel  for 
simulating  turbulence  transport  properties  in 
recirculating  flows  have  been  assessed  by 
comparing  continuum  model  predictions  of 
the  turbulent,  backmixed  flowfield  in  an 
opposed-jet  combustor  to  experimental  mea- 
surements. 

Both  eddy  viscosity  models  qualitatively 
describe  the  bulk  hydrodynamic  flowfield,  but 


the  detailed  flow  structure  is  inadequately 
represented.  Mass  transport  estimates  based 
on  the  algebraic  viscosity  model  conform  fa- 
vorably to  experiment.  The  two-equation 
turbulence  submodel  yields  consistently  poor 
correlation.  The  deficiency  of  the  two-equation 
eddy  viscosity  model  predictions  is  attributed 
to  the  boundary  condition  specification  for  the 
turbulence  energy  dissipation  rate.  It  is  shown 
that  the  turbulence  energy  dissipation  rate 
adjacent  to  critical  solid  walls  strongly  influ- 
ences the  overall  mixing  characteristics  of  the 
two-equation  submodel. 

The  results  indicate  that  the  algebraic  eddy 
viscosity  model  provides  cost-effective  predic- 
tions of  the  general  fluid  flow  patterns  and 
mass  transport  trends  in  turbulent,  confined 
flows  exhibiting  strong  recirculation.  The 
two-equation  eddy  viscosity  submodel  pro- 
vides better  resolution  of  small-scale  turbu- 
lence processes  but  requires  careful  testing  to 
ensure  realistic  predictions. 

It  is  also  shown  that  considerable  testing 
of  continuum  flow  models  in  isothermal  flow 
is  required  before  proceeding  to  the  compli- 
cating conditions  of  combustion.  The  hot  flow 
calculations  accentuate  the  inadequate  trans- 
[>ort  characteristics  identified  in  the  isothermal 
flow  analysis.  Attempts  to  refine  the  chemistry 
submodel  or  to  quantify  the  interaction  of  fluid 
motions  on  chemical  reaction  rates  are  depen- 
dent on  the  correct  turbulence  submodel  for- 
mulation. 

Nomenclature 

h stagnation  enthalpy 

k turbulence  kinetic  energy 

m mass  flowrate 

p pressure 

r radial  coordinate 

combustion  chamber  radius 
17,  V mean  velocity  components 

li ',0 ' fluctuating  velocity  components 

0 average  velocity 

X axial  coordinate 

Y mass  fraction 

r exchange  coefficient 

e turbulence  energy  dissipation  rate 

|x  viscosity 

p density 

CT  Prandtl  or  Schmidt  number 

<|>  dependent  variable 

<I>  equivalence  ratio 

Subscripts 

1 species 

in  inlet 
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j jet  stream 

m main  stream 

t turbulent 
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COMMENTS 


D.  B.  Spalding,  Imperial  College,  England.  I • 

would  like  to  mention  that  the  algebraic  viscosity 

model  is,  I think,  one  invented  by  W.  M.  Pun  and  GaryP.  Schneyer,  Science  Applications,  Inc.,  USA. 
myself  for  an  early  demonstration  of  numerical  I would  like  tof  ommend  the  authors  on  their  care 
combustion  models.  It  was  never  expected  to  have  wisdom  on  proceeding  first  with  isothermal 

wide  use.  That  it  has  survived  so  long  is  due  more  cases  and  only  then  proceeding  to  combustion.  I 
to  its  attractive  simplicity  than  to  its  validity  or  secure  have  a ccftnmeril  and  a question.  The  comment  is 
foundation.  that  for  many  “steady”  combustion  configurations, 

experimerfts  show  large-scale  unsteady  (often  peri- 
odic) motions.  For  your  configuration  I would  expect 
Authors'  Reply.  We  anticipated  that  the  algebraic  K-H  waves  to  appear  around  the  jet  causing 

viscosity  model  would  not  survive  the  challenge  asymmetries  if  nothing  else.  These  configurations 

to  experiment.  The  two  models  were  chosen  with  can  only  be  properly  modelled  by  transient  computer 

the  expectation  that  the  two-equation  model  would  codes.  Thus,  your  lack  of  agreement  may  be  affected 

neatly  resolve  the  deficiencies  produced  by  the  by  non-turbulent,  unsteady  phenomena  in  your  ex- 

algebraic.  That  this  did  not  occur  in  the  present  periments  so  that  your  turbulent  model  may  be  better 

experiment  should  cause  us  all  to  pause  and  ponder  than  you  think  (but  maybe  not).  ’ 

the  role  simplicity  should  play  in  our  methodology.  The  question  is  whether  you  measured  different 
Perhaps  additional  confrontation  with  experiment  turbulent  intensities  in  the  burning  vs.  isothermal 

shall  force  the  retirement  of  the  algebraic  model.  cases  that  were  otherwise  identical?  If  so,  how  did 

Perhaps  not.  they  differ? 
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I Authors' Reply.  At  present  there  is  no  experimental 

j evidence  of  large-scale,  transient  motions  in  the 
. opposed-jet  combustor  flowfield.  A laser  velocimeter 

' system  is  being  incorporated  into  the  OJC  test 
I facility  to  obtain  isothermal  and  hot  flow  turbulence 

j intensity  measurements  that  have  heretofore  been 

I unavailable.  This  information  will  broaden  the  data 

f base  for  evaluating  the  turbulence  submodels  and 

I will  also  facilitate  exploring  the  interaction  between 

[ fluid  motions  and  chemical  reaction  rates  in  turbu- 

[ lent,  backmixed  combustion. 


minimize  the  total  throughput  of  CO.  For  a maxi- 
mum CO  flowrate  of  Ig/sec,  concentrations  at  the 
exit  plane  are  less  than  those  measured  in  the 
products  of  combustion  for  the  reacting  case.  The 
venting  system,  designed  for  the  reacting  case,  is 
adequate  to  protect  experimenters  from  exposure  to 
CO. 

Tracer  studies  were  not  conducted  for  the  case 
of  hot  flow  in  the  present  study. 

• 


L 


R.  A.  Duerr,  NASA-Lewis  Research  Center,  USA. 
Carbon  monoxide  has  been  previously  used  as  a 
flow  tracer  in  a high  temperature  oxidizing  stream 
to  study  coaxial  jet  turbulence  where  the  CO  was 
converted  to  COj.  In  order  to  achieve  reasonably 
accurate  contour.;,  relatively  large  amounts  of  CO 
should  be  used  to  ensure  repeatability.  Assuming 
that  the  same  open  duct  apparatus  was  used  for 
cold  flow  as  well  as  hot  flow,  what  was  the  flow 
rate  of  CO  and  what  safety  precautions  were  taken? 
The  use  of  CO  as  a cold  flow  tracer  can  obviously 
be  rather  harmful  to  the  experimenters  over  an 
extended  test  period.  Other  gases  may  be  preferable 
from  the  standpoint  of  safety,  such  as  carbon  dioxide 
which  would  basically  involve  only  a change  in 
sample  cells  in  the  NDIR  analyzer  used  in  the  test. 

Authors’  Reply.  One  hundred  percent  (100%)  car- 
bon monoxide  (CO)  is  introduced  as  a tracer  species 
through  the  jet  for  the  isothermal  conditions.  Carbon 
dioxide  (COj)  is  used  as  a tracer  in  the  approach 
stream  to  allow  simultaneous  measurement  through- 
out the  flowfield  of  CO  (from  the  jet)  and  CO2 
(from  the  main  stream). 

CO  is  purposefully  selected  for  the  jet  tracer  to 


John  B.  Fenn,  Yale  University,  USA.  The  authors 
are  to  be  complimented  for  their  intrepid  attempt 
to  bridge  the  gap  between  measurements  on  and 
analysis  of  highly  idealized  combustion  models  and 
more  realistic  practical  burners.  I find  particular 
satisfaction  in  their  evidence  that  eddy  transport 
plays  a key  and  definable  role  in  the  behavior  of 
ducted  burners  and  provides  a basis  for  charac- 
terizing their  performance.  During  the  previous 
gathering  of  this  clan  at  MIT  my  then  colleague 
Emory  Wilkerson  and  I presented  some  results  on 
the  behavior  of  ducted  ramjet  burners.*  They  seemed 
clearly  to  show  that  measured  values  of  eddy  dif- 
fusivity  in  conjunction  with  a homogeneous  reaction 
model  for  the  combustion  process  were  remarkably 
able  to  describe  the  dependence  of  combustion 
efficiency  on  pilot  heat  input  and  burner  geometry. 
The  present  extension  of  those  then  primitive  efforts 
is  a progressive  stride  toward  the  goal  of  rational 
burner  design. 
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20  abstract  ^Conf^nuft  on  reverse  aide  It  neceaaary  and  Identify  by  block  nttmber) 

The  turbulent  flowfield  in  a backmixed  combustor  is  modeled  analytically  by 
simultaneously  solving  the  governing  partial  differential  conservation  equations. 
An  algebraic  and  a two-equation  eddy  viscosity  model  are  employed  in  the  numerica 
simulation  to  account  for  the  turbulence  transport  processes.  Predicted  distri- 
butions of  isothermal  flow  properties  are  systematically  compared  to  experimental 
obtained  data  to  appraise  the  eddy  viscosity  models.  In  addition,  predicted 
distributions  of  reacting  flow  properties  are  presen^Sd  to  illustrate  the  applica 
bllity  of  current  numerical  methods  to  the  predictio^  otSgont inuous  combustion.  . 
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The  <turbulent  momentum  and  mass  transport  properties  are  evaluated  in  isothermal 
flow  for  a range  of  mixing  conditions.  Botli  eddy  viscosity  models  qualitatively 
describe  the  system  hydrodynamics,  but  the  detailed  flow  structure  is  inadequate! 
represented.  The  mass  transport  prediction  from  the  algebraic  viscosity  mode 
agree  favorably  with  experiment.  The  inferior  performance  of  the  two-equation 
viscosity  model  is  remedied  by  refining  the  boundary  condition  specification  for 
the  turbulence  energy  dissipation  rate.  It  is  shown  that  the  turbulence  energy 
dissipation  rate  adjacent  to  critical  solid  walls  strongly  Influences  the  over- 
all mixing  characteristics  of  the  two-equation  model.  The  isothermal  flow 
results  indicate  that  the  algebraic  eddy  viscosity  model  provides  cost-effective 
predictions  of  the  general  fluid  flow  patterns  and  mass  transport  trends  in 
confined  flows  exhibiting  strong  recirculation.  The  two-equation  eddy  viscosity 
model  provides  better  resolution  of  the  small-scale  turbulence  processes  but 
requires  careful  testing  to  ensure  realistic  predictions.  The  hot  flow  calcula- 
tions accentuate  the  inadequate  transport  characteristics  identified  in  the 
isothermal  flow  analysis  and  varify  that  considerable  testing  of  the  numerical 
model  is  required  before  proceeding  to  the  complicating  conditions  of  combustion. 
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